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Abstract 

Stationary equilibria of point vortices with arbitrary choice of cir- 
culations in a background flow are studied. Differential equations sat- 
isfied by generating polynomials of vortex configurations are derived. 
It is shown that these equations can be reduced to a single one. It 
is found that polynomials that are Wronskians of classical orthogonal 
polynomials solve the latter equation. As a consequence vortex equi- 
libria at a certain choice of background flows can be described with 
the help of Wronskians of classical orthogonal polynomials. 

1 Introduction 

The model of point vortices describing motion of two-dimensional incom- 
pressible fluid is one of the most elegant models of fluid dynamics. The 
motion of M point vortices with circulations (or strengths) T\, ... , Tm at 
positions z±, . . . , Zm in zero background flow is governed by the Helmholtz's 
equations 

d -^ = — Y'-^—, k = l,...,M, 

j=i j 

where the prime means that the case j = k is excluded and the symbol * 
stands for complex conjugation. This system induces complicated dynam- 
ics and is not integrable in the case M > 3. Nevertheless some particular 
"motions" including stationary, translating, and rotating equilibria are suc- 
cessfully studied PHI]. A convenient approach applicable to such types of 
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motion is the so-called "polynomial method" [2J. According to this method 
polynomials with roots at vortex positions are introduced. This approach 
provides quite unexpected connection between dynamics of point vortices and 
the theory of classical and nonlinear special polynomials [2|IT2 | ll3 pi5l4T8] . For 
example, the generating polynomial of M identical point vortices on a line 
in rotating relative equilibrium is essentially the Mth Hermite polynomial. 
In 1964 Tkachenko obtained a differential equation satisfied by generating 
polynomials of stationary vortex arrangements with equal in absolute value 
circulations [19]. Now this equation is known as the Tkacheno equation. 
Later a generalization of the Tkachenko equation to the uniformly trans- 
lating case was derived (for details see [USE]). Again the latter equation 
describes translating equilibria of point vortices with equal in absolute value 
circulations. In recent work [TJ] it is shown that equations for generating 
polynomials of vortex arrangements with arbitrary choice of circulations can 
be reduced to the Tkachenko equation in the completely stationary case 
and to the generalization of the the Tkachenko equation equation in the uni- 
formly translating case. It turns out that the Adler - Moser polynomials [15J, 
a famous sequence of polynomials, which generates polynomial solutions to 
the Tkachenko equation, provide not unique polynomial solutions of the lat- 
ter equation. Some examples of alternative polynomial solutions are given 
in [H]. 

In this article we study stationary equilibria of multivortex systems in 
a background flow. We derive differential equations satisfied by generating 
polynomials of vortices and reduce these equation to a simple form. We ap- 
ply an approach suggested in rfj] . Our aim is to show that vortex equilibria 
at a certain choice of background flows can be described with the help of 
polynomials that are Wronskians of classical orthogonal polynomials. These 
results provide additional link between the vortex theory and the the theory 
of classical orthogonal polynomials. We use the technic of Darboux transfor- 
mations, for more details see p~5lf20]. l2Tj . 

This article is organized as follows. In section |2j we construct sequences of 
Darboux transformations for a second order linear differential equation and 
study in details the cases corresponding to the families of classical orthog- 
onal polynomials. In section [3] we consider stationary equilibrium of point 
vortices in a background flow and derive differential equations satisfied by 
generating polynomials of vortex configurations. We give background flows 
in explicit form for such cases that involve Wronskians of classical orthogonal 
polynomials. In addition in section [3J we construct reductions of the differen- 
tial equations satisfied by generating polynomials to a simple one. In section 
H] we study polynomial solutions of the latter equation and present several 
explicit examples. 
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2 Darboux transformations for classical or- 
thogonal polynomials 

The systems of classical orthogonal polynomials can be constructed as poly- 
nomial solutions of the following second order linear differential equation 

(j{z)i) zz + t(z)^j z + \i) = 0. (1) 

In this expression a(z) is a polynomial of the degree at most two, r(z) is a 
polynomial of the degree at most one, and A is a real constant. First of all 
we shall consider a more general second order equation 

v(z)ip zz + T(z)ip z + {u(z) + A}-0 = (2) 

with a(z), t(z), u(z) being sufficiently smooth functions. Let us denote a 
nontrivial solution of equation (J2J) with the parameter A = Ai as ipi. 

Theorem 2.1. For any solution ip of equation (TS]) the Darboux transfor- 
mation 

4> = i> z - ^V>, (3) 
gives a solution of the following equation 

o-(z)^ zz + n(z)4 + {ui(z) + = 0, (4) 

where the functions Ti(z), u%(z) are given by 

n(z) =t(z) +a z (z), Ux{z) = u(z) + r z {z) + a z (z){lnip 1 } z + 

+2a(z){\mp 1 } zz . 

Proof. We introduce the second order operators L and L\ according to 
the rules 

L = f o-(z)D 2 +t(z)D+u(z), L 1 = f a(z)D 2 +T l (z)D+u l (z), D d = ^- (6) 

and rewrite equations 02]), dl]) in the following way 

(L + = 0, (Li + \)4> = 0. (7) 
Substituting the transformation ([3]), which we rewrite in the form 



into the second equation in (j7|) and using the first equation in (j7|), we obtain 
the relation 



(£iA/>i - A^Ljiff = 0. (9) 

Expression (jUj) is a first order polynomial in ip zz , ip z , ip. Setting to zero coeffi- 
cients at ipzz, ipz yield relations (0) for the functions Tx(z), u\(z) accordingly 
Using these relations and the equation for the function ipi, we see that the 
coefficient at ip in expression (j^J) vanishes. This completes the proof. 

Further we note that the Darboux transformation can be iterated. As 
a result we obtain the following theorem. 

Theorem 2.2. Letipi, . . ., ip k be nontrivial solutions of equation (j2J) with 
pairwise different values of the parameter A: Ai, . . ., X k . Then the Darboux 
transformation 

f = wbh^ks (10) 

5W es a solution of the equation 

°{z)i) zz + r k {z)i) z + {u fc (z) + A}^ = (11) 

def 

whenever the function ip solves equation (j2J). In this expressions W = 
W[ipi, . . . , ipk] is the Wronskian of the functions ipi, . . ., ip k and 

T k (z) = t(z) + ka z (z), 

k(k — 1) 

u k (z) = u(z) + kr z (z) + ^— — L a zz (z) + a z (z){ln W} z + (12) 
+2a(z){lnW} zz . 

Proof. Applying k times the Darboux transformations to a solution ip 
of equation (jSJ), we see that the function 

?P = M k ij, M k = f A fk ...A f2 A fl , h = ^ (13) 
satisfies equation ffTTl) with r k (z), u k (z) given by 
T k (z) = t(z) + kcr z (z), 

k(k — 1) 

u k (z) = u(z) + kr z (z) + ^— — L o zz {z) + a z (z){ln h . . . f k } z + (14) 
+2a(z){lnf 1 ...f k } zz . 
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We note that at the n th step, n — 1, . . ., k we use the Darboux transformation 
build on the basis of the function f n . Thus we obtain the set of the functions 

/l> ■ ■ •> fk 

fi = fa, h = A h ip 2 , ■ ■ ■ , fn = A fn _, . . . A h A fl ip n , f k = 

= A f k -i ■■■ A h A h i>h 

From our construction it follows that the functions ipi, . . ., ip n provide n 
linearly independent solutions of the n th order linear differential equation 

M n w = 0, M n = f A fn ...A f2 A fl . (16) 

Further we rewrite this equation in explicit form 

D n w+m n _ 1 D n - 1 w + . . .+m w = 0, M n = D n +m n _ 1 D n - l + . . .+m . (17) 

Substituting the functions fa, . . ., fa into equation (ITT]) , we get the system 
of linear algebraic equations with respect to the coefficients mo, . . ., m„_i 

m fa + m x Dfa + ... + 77in_iD n -Vi = ~D n fa l = l,...,n. (18) 
Solving this system with the help of the Kramer's rule yields 

Wi 

m-i = -TTT7- -7 / = l,...,n, (19) 

W[fa,...,fa] 

where the determinant W\ is obtained by replacing the Z th column of the 
Wronskian W[fa, . . . , fa] by the column of the functions D n ipi, . . ., D n ip n . 
Consequently, we see that the operator M n acts on any function w as 

= Wfr 4, *] 
W[fa,...,fa] 



Now we return to Darboux transformation (IT3]) . It follows from relation f l20|) 
that transformation ffl3l) is given by fflQj) . Expressions (ITT]) are obtained by 
induction applying k times formulae (J5J). It remains to find the product f\ 
. . . fk- Using expressions ( 1T5|) . (TTB|> . (12Dj). we obtain 

fi- ■■//, = * . . . ^rr^i = w ^ ■■■■*]■ < 21 > 

V>i VK [-01, • • • , Vv-ij 



Substituting relation (T2T]) into expressions (JHJ), we see that the function ([10 
is a solution of equation (ITT]) provided that conditions ( TT2]) hold. 
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Now let us apply theorem 2.2 to the equation for classical orthogonal 
polynomials (see (□])). Suppose p^, . . ., p ik , Pi k+1 are classical orthogonal 
polynomials satisfying equation (JTJ) with pairwise different values of the pa- 
rameter A: Aj x . . ., Ai fc , Aj fe+1 . Then the polynomials 



P( z ) = W\p h , . . . ,Pi h ,p ih+1 ], Q(z) = W\p h: 
solve the bilinear equation 



a {P ZZ Q - 2P Z Q Z + PQ ZZ } + 



k — - I rr. 



{P Z Q - PQ Z } 



(22) 



(7, 



+-£{P.Q + PQz} + 



k(k-i] 



Ozz + kr z + A 



(23) 



PQ = 0, A = A, fc+1 . 



This result is obtained substituting ip = P/Q, u(z) = 0, and expressions (TT2"j) 
with W = Q(z) into equation (TTTT) . where we set A = X ik r Further finding 
the highest powers of Wronskians, we obtain the degrees of the polynomials 

P(z),Q(z) 



degP = ^ l -^^ 1 , de g Q = J2*i- 



k(k-l) 



(24) 



i=i 



i=\ 



where i\ is the degree of the polynomial Pi r In the next section our aim is 
to show that equation (1231) can be obtained in the framework of the vortex 
theory. 



3 Point vortices in a background flow 

The motion of a multivortex system in a background flow w(z) is described 
by the Helmholtz's equations with additional term 

dz l - 1 y> T _i | w *^\ k = l,...,M. (25) 

d t 2iri ^— ' zi, — Za 2iri ' 
j=i j 

In what follows we take the background flow in the form 

U(z k ) T k R z (z k ) 
W W = 2RW) + TflW- (26) 

where R(z) is a polynomial of degree at most two and U (z) is a polynomial 
of degree at most one. We exclude the case when R{zk Q ) = unless z = z ko is 
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a simple root of the polynomial R(z) and U(z ko ) + T ko R z (zk )/2 = 0. Let us 
subdivide the vortices into groups according to the values of their circulations 
Tj, j = 1, . . . , N. In other words, we suppose that there are N different 
values of circulations in the arrangements we study. By , . . . , a| we 
denote the positions of the vortices with circulation Tj. Therefore, we have 
M = li + . . . + l N . Let us introduce the polynomials 

h 

P 3 (z) = \\(z - a!f\ j = l,...,N. (27) 
i=i 

with roots at the vortex positions. We shall consider the stationary case, 
this yields dz^/dt = 0. From equations (125]) we find the system of algebraic 
relations 

V V' Ez - ( a ) i Rg ( Q ) = n 

a-af ZR{o) 4 it! (a) ' (28) 

a = a S o) > *o = 1, • • • , lj , jo = l,...,N, 

where the case (jo,io) = in summation is excluded. Using properties of 
the logarithmic derivative, we obtain the following equalities 



,Oo) 

lating the limit z — > a in the expression for Pj Q ZZ yields 



Now let z tend to one of the roots a = a\ J o ° of the polynomial Pj (z). Calcu- 



In 



P J(hZZ (a) = 2P io »^' fry- (30) 

i=i a - a i 

Using equalities f[28]) . f )29|) . we get the conditions 

r Pj ,zz(a) _ 9 PjM U (a) r j0 R z (a) 

i0 P f0 ,(a) ^ j P(a) P(a) 2 P (a) ' 1 J 



which are valid for any root a = af° of the polynomial Pj (z). Further we 
see that the polynomial 

N { N P P P 

1 ' '' ' " 3 (32) 



being of degree M + max{deg R, deg U + 1} — 2 possesses the same roots as 
the polynomial Pi (2) x ... x Pj^(z). Thus we find the equation 



j=l 3 k<j 

n , r 2 \ 

.7=1 ^ ' 



p V / / j K J r> p 

■ "j r~ P fe "7 

(33) 

AT , -n 2 n p V7 

where the constant C is obtained setting to zero the coefficient at z° in the 
Laurent expansion of the left-hand side in the relation (1331 : 

c = ^P_f^ VjdBgP \ +Uz(z) jrr j degP j . (34) 
\j=i / j=i 



Further following an approach suggested in the article fTJ] , we introduce new 
functions P(z), Q(z) according to the rules 



N r ,( r ,+ 1 ) _ JL r j( r j~ 1 ) 



P(z)=l[P J 2 (z), Q(z) = l[P J 2 (z). (35) 
i=i j=i 



Using the equalities 



|>{%)<5m} = £i1 

j=i \ 



" 'p. Pr 

J ],ZZ ± J,Z 



p. P2 / ' 



'JSflUfr,^, <36) 

we get the differential equation for the functions P(z) and Q(z) of the form 
rf 2 , r~, si . ( d , i P(z))Y „, , d , fPM 



Note that this equation can be rewritten in the bilinear form as 

{PzzQ - 2p z q z + pg 22 } + u(z) [p z q - pg,} + 

8 



(37) 



(3? 



If P, Q in (|38p are polynomials, then for the constant C we obtain 

C = (deg P - deg Q) 2 + U g (z) (deg P - deg Q) . (39) 

Thus we see that stationary equilibria of point vortices with circulations r 1; 
. . ., Tjv in the background flow given by (1261) can be described with the help 
of equation (1381) . Suppose we have found a solution of equation (1381) in the 
form (|35|) . then a vortex with circulation Tj is situated at the point z = zq 
whenever the function P(z) has a "root" of "multiplicity" Tj(Tj + 1)/2 at the 
point z = z and the function Q(z) has a "root" of "multiplicity" Fj(Tj — l)/2 
at the point 2 = zq- The circulation is calculated as the difference of the 
corresponding "multiplicities". We say that the function f(z) has a "root" 
of "multiplicity" k at the point z = z if /(z) = (2 — z ) k ip(z) with ^(z) 
being an analytic function (possibly multivalued) in a neighborhood of Zo 
and ip(z ) 7^ 0. Suppose /c G N, then the point z = z is a root of the 
function /(z) in the usual cense. 

Now we shall apply our results to the following multivortex situation. 
Suppose /+ vortices with circulations nT are situated at positions z = a[ n \ 
i = 1, . . ., l£, n = 1, . . ., Ni and /~ vortices with circulations —mT are 
situated at positions z = b^, k = 1, . . ., m = 1, . . ., N 2 . Further we 
consider the polynomials 

1+ 

P n (z) = l[(z-a ( [ l) ), n=l,...,iV 1 

l=1 _ (40) 

Q m (z) = H(z-b^), m = l,...,N 2 



fc=i 



with no common and multiple roots. The amount of vortices in such an 
arrangement is equal to M = If + . . . + it + V{ + . . . + l^. We see that 
stationary equilibria of this multivortex system in the background flow 
can be described in terms of the functions 



1 rn(rn+i) \ rm(rm-l) 



n=l m=l / 4 -,n 

JVi V9 v ' 

1 rn(rn-i) * rm(rm+i) 



P(z) = n^n 2 WHQm 2 (Z) 

n=l m=l 

n=l m=l 

which satisfy equation ( I38p . Along with this we can consider the polynomials 
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1 n(n + l) ^ 2 m(m-l) _ ^\ n(n-l) ^ 2 m(m+l) 



(42) 



(43) 



of the form 

p{z) x\r—( ) n qt^(z), q{z) = i[Pn * (z) n q 

n=l m=l n=l m=l 

and conclude that these polynomials satisfy the equation 

R(z) [P ZZ Q - 2P Z Q Z + PQ 22 } + ^ {P Z Q - PQ Z } + 

From expression (|39|) it follows that 

C: = (degP - degQ) 2 + ^ (deg P - de g g) . (44) 

Now recalling the results of section [2J we obtain that equation (1381) coin- 
cides with equation f[2"3"l) provided that 

R(z) = a(z), U(z) = t(z) + (k - <t*(z) (45) 

and -P(-z) = P(^), = By direct calculations we verify that the 

coefficient at PQ in (1381) coincides with the coefficient at PQ in (1231) whenever 
relation (|45p is valid. Consequently, the polynomials 

P(z) = W\p h , . . . ,p ik ,p ik+1 }, Q(z) = W\p h ,. . ., Pik ] (46) 

where p ix , ... , Pi k+1 are pairwise different classical orthogonal polynomials, 
solve equation (j58|) provided that relations fl4"5]) hold. 

4 Hermite, Laguerre, and Jacobi polynomials 
in the framework of the vortex theory 



In this section we study polynomial solutions of equation (1381) and give several 
explicit examples. We begin with the theorem. 

Theorem 4.1. Suppose a pair of polynomials P{z), Q(z) satisfy equa- 
tion (1351) and z = z is such a point in the complex plane that R(z ) ^ 0; 
then the following statements are valid: 
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1. If the point zq is a multiple root of one of the polynomials P(z), Q(z), 
then it is also a root of another polynomial and the multiplicities of the 
root zq for the polynomials P{z), Q(z) are two successive triangular 
numbers. 

2. If the point Zq is a common root of the polynomials P(z), Q(z), then it 
is a multiple root of at least one of them and the multiplicities of the 
root zq for the polynomials P{z), Q(z) are two successive triangular 
numbers. 

Proof. In order to prove statements 1, 2 of the theorem we substitute 
the expressions 

P(z) = {z- z ) r h{z), Q(z) = {z- z ) s g(z), r,seNU {0}, (47) 

where h(z), g(z) are polynomials such that h(zo) ^ 0, g(zo) ^ into equation 
(I38p and find the Tailor series in a neighborhood of the point Zq of the 
resulting relation. Setting to zero the coefficient at (z — z ) r+s ~ 2 , we obtain 
the algebraic equation 

(r - s) 2 = r + s. (48) 
Thus we get the system 

r-s = l, r + s = l 2 , I el. (49) 
Solving this system yields 

r = !G±i>, . = ^, leZ. (50) 

Analyzing expressions fl50l) . we make sure that statements 1, 2 are valid. 

As a consequence of theorem 4.1 we see that any polynomial solution 
P(z), Q(z) of equation fl38|) such that the polynomials P(z), Q(z) do not 
have common roots with the polynomial R(z) describes equilibrium of point 
vortices in a background flow. Indeed, from theorem 4.1 it follows that the 
polynomial solution in question can be always presented in the form (1551) . 

First of all let us study the case when one of the polynomials (P(z) or 
Q{z)) in (!38|) is a constant. Setting Q(z) = 1 in ([58"]) . we obtain the equation 
for the polynomial P(z) 



R{z)P z 



^ + U(z) 



P Z -CP = 0, (51) 
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(b) LR^(z) 



(c) LRf\z) (d) LRf\z) 

Figure 1: Roots of the polynomials {LRj? (z)}. 



where the constant C is given by with degQ = 0. The equation ( 15Tp 
coincides with the equation for classical orthogonal polynomials (see (pQ)) 
provided that the polynomials R(z), U(z) are taken in the form (|45|) with 
k = 0. It is well known that classical orthogonal polynomials do not have 
multiple roots and do not have common roots with the polynomial R(z). 
Thus it follows from expression (j35p that the roots of any classical orthogonal 
polynomial give equilibrium positions of M — degP point vortices with 
circulation 1 in background flow f l26|) with R(z), U(z) taken as ( H5j) under 
the condition k = 0. Analogously setting P(z) = 1 in fl38|) . we obtain the 
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Table 1: Polynomials {LRf'jz)}. 

LR [ a) (z) = l 
LR[ a \z) = z-a-1 

LR ( 2 a \z) = z 3 - 3{a + 2)z 2 + 3(a + l)(a + 3)z - (a + l)(a + 2)(a + 3) 
LRf\z) = z % - 6(a + 3)z b + 15(a + 2) (a + 4)z 4 - 10(a + 3) (2a 2 + 12a + 13)^ 3 
+15(a + l)(a + 3) 2 (a + h)z 2 - 6(a + l)(a + 2)(a + 3)(a + 4)(a + 5)z 
+ (a + l)(a + 2)(a + 3) 2 (a + 4) (a + 5) 



equation for the polynomial Q(z 

'RAz 



R{z)Q zz + 



U(z) 



Q Z -CQ = 0, (52) 



where the constant C is given by (|3"U|) with degP = 0. This equation is 
exactly the equation for classical orthogonal polynomials (see ([TJ) if we take 
the polynomials R(z), U(z) in the form 

R(z) = a(z), U(z) = 2&-t(z). (53) 

Consequently from expression (|35p it follows that the roots of any classical 
orthogonal polynomial give equilibrium positions of M — deg Q point vortices 
with circulation —1 in background flow (|26|) with the polynomials U(z), R{z) 
given by (1531 . 

Further recalling the results of section [21 we consider polynomials that 
solve equation (1381) and are Wronskians of classical orthogonal polynomi- 
als. The sequence of Hermite polynomials satisfy the following second order 
differential equation 

ij zz - 2zt/j z + 2nij = 0, tfj = H n {z). (54) 

Thus, we have a(z) = 1, t(z) = —2z, X n = 2n. The Hermite polynomials 
are orthogonal with respect to the weight function q{z) = exp(— z 2 ) on the 
real line (—00, 00). Using expressions (1381) . (145|) . we see that the equation 

P ZZ Q - 2P Z Q Z + PQ ZZ - 2z {P Z Q - Pg 2 } + 2 (degP - degQ) PQ = 0. 

(55) 

possesses polynomial solutions given by formula (j4"6|) with p ik = H ik (z). 
Originally this result appeared in [21]. It is known that polynomials that 
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-0.5 



0.5 

a Rez 



(b) UG 7 (z) 



-0.5 



(c) UGs(z) (d) UG 9 (z) 

Figure 2: Roots of the polynomials {UGk{z)}. 



are Wronskians of the Hermite polynomials also arise in the theory of the 
Painleve equations and their higher order analogues [12, 22 J. 

The sequence of Laguerre polynomials L„ (z) can be generated with the 
help of the following ordinary differential equation 

zip zz + (a + 1 - z)i) z + nip = 0, V = L ( n a) (z), a > -1. (56) 

Note that sometimes these polynomials are called the generalized Laguerre 
polynomials whenever a / 0, a > —1. In this case we get a(z) = z, t(z) = 
a + 1 — z, \ n = n. The Laguerre polynomials are orthogonal with respect to 
the weight function q(z) = z a exp(-z) on the real interval [0, <x>). By means 
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Table 2: Polynomials {UGk(z)}. 



UG t {z) =z z -\z 

UG 2 {z)=z*-lz* + ^z 2 + ± 1 

UG 3 (z) = z* - \z 7 + lz 5 + l 4 z 3 - §.z 

jjn (A — ? i2 _ i r u ± i ? s ^ I ? e _ i a j , _ls_ 

UUrzt^j ^ 2 4, T T 25 g^ T 25 g^ T 40g6 

Tin — r 15 -I- liir 11 -I- Ii>r 9 — 165 -,7 i 147 r 5 _ _35_ r 3 _ 105 _ 

^<^5l^J — *> 16 "^16^ 256^ "^512^ 4096^ 8192^ 

Tin (~\ _ -.18 I 3 16 I 15 .,14 I 175 r 12 231 10 _i_ 231_ 8 49 -6 

^ "I" 4^ " r 8 41 " r 64 ^ 256^ " r 256 x; 2048^ 

735. 4 , 735 2 , 49 



8192 ~ 1 32768 1 65536 



of expressions (155)1 . we obtain that the equation 

z |P 22 Q - 2P Z Q Z + PQ«} + (a + k + ^-zj [P Z Q - PQ Z 
+1 {PzQ + PQ Z ) + (cleg P - degQ) PQ = 0. 



Li#>(z) 



(57) 



has polynomial solutions given by formula (l46p with p ife = L^(z). Let us 
consider an example. Two neighbor polynomials from the sequence 



(5? 



i.e. P(z) = LR^^z), Q(z) = LR^*\z) solve equation (jSTj) . This statement 

for the pair P(z) = LR^[\z), Q(z) = LRq(z) follows from equation fl51"|) 
and remarks after it. Without loss of generality, we choose the constant Sk 
in such a way that the resulting polynomial is monic. First few polynomials 
from the sequence (j5"8]> are given in tabled) Using relation f[2"4l) . we find the 
degree of the polynomial LR^\z) 

degL4^(,) = ^±ll. (59) 

Interestingly, that roots of the polynomials LR k a \z) form highly regular 
structures in the complex plane. Examples of plots are given in figure (T) 
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Figure 3: Roots of the polynomials {LGrfc(z)}. 



Classical polynomials orthogonal on a finite interval, which is usually 
taken as the real interval [—1,1], are the Jacobi polynomials Pn*'^\z) and 
their partial cases including the Chebyshev polynomials of the first kind T n (z) 
and of the second kind U n (z), the Gegenbauer polynomials Cn (z), and the 
Legendre polynomials P n (z). The Jacobi polynomials satisfy the following 
second order equation 

(1 - z 2 )i) zz + [(3-a-(a + /3 + 2)z) ip z + n(n + a + (3 + l)tp = 

^ = Pi a ^(z), a >-l, (3>-l. (60) 

Consequently, we have a(z) = 1 — z 2 , t(z) =(3 — a — (a + (3 + 2)z, X n = 
n(n + a + (3 + 1). Note that the Jacobi polynomials are orthogonal with 
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Table 3: Polynomials {LG k (z)}. 



LGi 




= z 




LG 2 




_ ~4 2 2 1 

* 7^ 35 




LG 3 




_ r 9 76 r 7, 114 5 1 4 3 1 _ 
~~ * 91^"^ 715^ "r" 143^ 143^ 




LG 4 




_ 16 2840 14 , 45964 12 856 10 
* 1729 ' 46189 4199 
140 ^4 , 1400 ^2 , 35 

31603 ~ 1 6605027 ~ 1 6605027 


1970 ^8 1 16408 ^6 
46189" 1 600457 



respect to the weight function q(z) = (1 — z) a (l + z) 13 on the real interval 
[—1,1]. Using expressions (138]) . (HSJ), we find the differential equation 

(1 - z 2 ) \P ZZ Q - 2P Z Q Z + PQ 22 } + 
+ \p - a - (a + (3 + 2k + l)z) {p 2 Q - PQ 2 } - z {p 2 Q + PQ 2 } + (61) 
+ (degP - degQ) (degP - degQ + a + + 2k + l) PQ = 0, 

which possesses polynomial solutions in the form f )46]) with pi k = P^^\z). 

In the case a = (3 = 1/2 the Jacobi polynomials Pn a '^\z) become the Cheby- 
shev polynomials of the second kind U n (z). As an example we consider the 
following sequence of polynomials 

UG k (z) = lk W[U 3} U 4 ,..., U k+2 ], fceN. (62) 

Two neighbor polynomials from this sequence, in other words the polynomials 
P(z) = UG k+1 (z), Q(z) = UG k (z), satisfy equation with a = (3 = 1/2. 
In expression fl62l) we choose the constant ^ k in such a way that the corre- 
sponding polynomial is monic. With the help of relation (|24l) we calculate 
the degree of the polynomial UG k (z) 

degUG k (z) =3k. (63) 

First few polynomials UG k {z) in explicit form are given in table [2] Again 
the roots of the polynomials UG k (z) form highly regular structures in the 
complex plane. Several examples are plotted in figure [2J 

The Legendre polynomials P„(z) are special cases of the Jacobi polyno- 
mials and correspond to zero values of the parameters a, (3. Let us consider 
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the following sequence of polynomials 

LG k {z) = u k W[P lt P^..., Psfc-2], keN. (64) 

Two neighbor polynomials from this sequence P(z) = LGk+i(z), Q(z) = 
LGk(z), satisfy equation ( I6"T|) with a = (3 = 0. In expression (jM]) we take 
the constant Uk in such a way that the corresponding polynomial is monic. 
By means of relation (T24"j) we derive the degree of the polynomial LGk{z) 

deg LG k (z) = k 2 . (65) 

Several polynomials LGk{z) in explicit form are presented in table [31 The 
roots of the polynomials LGk(z) form highly regular structures in the complex 
plane. Several examples are given in figure [31 

Finally, we would like to mention that our investigation establishes an 
additional connection between the point vortex theory and the theory of 
classical orthogonal polynomials. We hope that our results will be useful in 
further studying of point vortex equilibria. 

5 Conclusion 

In this article we have studied the problem of finding stationary configura- 
tions of point vortices with generic choice of circulations in a background flow. 
We have found differential equations satisfied by generating polynomials of 
vortex arrangements and have shown that these equations can be reduced to 
a single one. We have studied the latter equation and have derived polyno- 
mial solutions expressed as Wronskians of classical orthogonal polynomials. 
This result was obtained with the help of an approach based on the tech- 
nique of Darboux transformation. In details we considered several examples 
involving Hermite, Laguerre, Chebyshev, and Legendre polynomials. 
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